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We use a cooling algorithm based on an improved action with scale invariant instanton solutions, which needs 
no monitoring or calibration and has a inherent cut off for dislocations. In an application to the SU{2) theory 
the method provides good susceptibility data and physical size distributions of instantons. 
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1. THE IMPROVED COOLING 

In we presented an improved cooling algo- 
rithm which fulfills requirements qualifying it as 
a method to eliminate UV noise and study topo- 
logical excitations on the lattice: 

- it smoothes out the short range fluctuations, 
including "dislocations" 

- it ensures stability of instantons (including their 
size) above a certain threshold po ~ 2.3a and 
needs thus no monitoring or calibration. 

We apply this method to SU(2) to study suscep- 
tibility, size distribution, distances and shape of 
instantons and we report here results concerning 
scaling and boundary condition effects. 

Since the Wilson action decreases with the in- 
stanton size, it leads to an unphysical abundance 
of small instantons in a Monte Carlo simulation. 
Under Wilson cooling even wide instantons shrink 
and decay. This can be corrected for by modify- 
ing the lattice action. Starting from |^ one can 
construct a one parameter set of actions with no 
0{a'^) and 0{a*) corrections using five fundamen- 
tal, planar loops of size m x n: 
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*Talk by Margarita Garcia Perez at LATTICE96. 



Here (m„ n,) = (1, 1), (2, 2), (1, 2), (1, 3), (3, 3) for 
i — 1, ... ,5 and 

ci = (19 - 55 C5)/9, C2 = (1 - 64 c5)/9 

C3 = (-64 + 640 C5)/45, C4 = l/5-2c5 (3) 
Our choice for an improved action is defined with 
the tuning cs = 1/20. The topological charge is 
improved the same way. We shall call this loop 
combination 5 Li: it will be exclusively used in the 
following. As shown in [|l| it ensures scale invari- 
ance for instantons above a threshold po ~ 2.3a. 
The cooling algorithm exactly minimizes the local 
action at each step and involves no further cali- 
bration or engineering. The size dependence of 
the action for the instantons is 9 - function like, 
giving Stf^ to better than 0.1% for p > po and a 
steep descent (slightly configuration dependent) 
in the interval (0.8 — 1.0)pa. Since the threshold 
is fixed in lattice units, it should not affect phys- 
ical sizes if a is small enough. For definiteness we 
use everywhere as radius in lattice units: 

P"* = Pteak = 6/ (TT^Qpeafe) (4) 

with Qpeak the topological charge density at the 
center of the instanton (using Speak instead gives 
no significant difference). Notice that the Q{5Li) 
charge operator produces reasonable density data 
already on still rough configurations and can 
be used to observe instanton-anti-instanton (lA) 
pairs (which annihilate later in the cooling since 
they are not minima of the action). 
The question of stability of instantons cannot be 
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completely disentangled from finite size effects. It 
is proven that Q = ±1 solutions do not exist on 
finite periodic volumes |^ . Practically this means 
that isolated instantons above ~ 1/ 4th of the lat- 
tice size are unstable. To avoid this affecting the 
physical distances one can use twisted boundary 
conditions or make at least one of the lattice sizes 
very large — both these conditions allow to sta- 
bilize wider instantons |^ . The results of were 
obtained using twisted b.c. with k — (1,1,1), 
m = (0,0,0) (t.b.c. in the following). 

The analysis undertaken here concerns the de- 
termination of physical quantities like topological 
susceptibility and size distribution in view of the 
two problems raised above, namely (A) threshold 
effects (for rough lattices po may already repre- 
sent a physically relevant distance) and (B) finite 
size effects (small lattices may not accommodate 
large, physical instantons). Correspondingly we 
analyze here the scaling behavior and the depen- 
dence on size and boundary conditions. 

2. RESULTS AND DISCUSSION 

We base our analysis on the following SU{2) 
Monte Carlo simulations (heat-bath, Wilson): 

(1) : 123xiVt, (3 = 2.4, (a ~ 0.12fm) (a): Nt = 12, 
t.b.c, (b): Nt = 12, p.b.c, and (c): Nt = 36, 
p.b.c; 

(2) : 12^ /? = 2.5, (a ~ 0.085fm), t.b.c; 

(3) : 24^, p = 2.6, (a ~ 0.06fm), t.b.c. 

Here (1) and (3) are new data from 200, re- 
spectively 84 configurations (100, respectively 200 
sweeps apart, after 20000 sweeps for thermaliza- 
tion), (2) are older data from 160 configura- 
tions (250 sweeps apart). 

Topological Susceptibility: Since the topological 
charge Q{5Li) stabilizes very fast to an integer 
within less than 1% (see Fig. 1), we expect devia- 
tions from the physical values in the susceptibility 
to be only of type (A) or (B) above. The results 
are presented in Fig. 2 and Table 1. The sus- 
ceptibility settles early in the cooling and seems 
to scale very well both with the cut-off (compare 
the (1) and (3) data) and with the volume (com- 
pare (1) and (2)), with small deviations show- 
ing up in (1) (beta=2.4) which will be argued to 
be mostly of type (A). Assuming the Witten - 



Veneziano formula to hold the result agrees ex- 
cellently with the phenomcnological expectation. 
Recent results for SU(2) in Q are about 15% 
higher than the ones we quote. 
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Figure 1. Evolution of action and charge of a 
typical configm'ation (t.b.c) under 5Li cooling. 
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Figure 2. Topological Susceptibility: t.b.c. (la), 
(2) and (3) (squares), and p.b.c. (lb) (triangles) 
and (Ic) (circles) - see text. The bars in (1) give 
an indication of the errors which are not included 
in the data to make the figure clearer. 
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Table 1: \/x (MeV) vs number of cooling sweeps. 



Instanton Size Distribution: The size of the in- 
stantons is calculated following the continuum 
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ansatz. We verified that a discretized form of this 
ansatz with periodicity corrections and general- 
ized to an elHpsoid fits very well the 5Li-densities 
(charge and action) of isolated instantons already 
after few cooling sweeps and that it can also han- 
dle overlapping structures. The fitted radius av- 
eraged over directions agrees quite well with ppeak 
Eq.(4). Wide instantons above 1/2 A^s show con- 
siderable anisotropy and for such configurations 
different definitions of p may lead to different val- 
ues. We have restricted to the Ppeak definition ex- 
cept for the old data in where the size was ex- 
tracted differently. Consequently the upper part 
of the size distribution for (2) [Q cannot be di- 
rectly compared with (1) and (3). 
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Figure 3. Size distributions. Crosses, squares, tri- 
angles, circles correspond respectively to lattices 
(3), (la), (lb), (Ic), the dotted vertical line is po- 



The new results are presented in Fig. 3. The 
comparison of the (1) and (3) data shows that 
at /3 = 2.4 the size distribution already begins 
to get out of the threshold region, indicating a 
reduced effect of the latter on physical instantons. 
The slight decay of the susceptibility observed in 
(1) can be completely understood as due to the 
disappearance of few, small but still physically 



relevant instantons close to the threshold. 

Finite size effects can be estimated by compar- 
ing (la,b,c). The results agree perfectly for low 
cooling sweeps but for (lb) (12^, p.b.c) a shift 
towards smaller sizes at 300 sweeps is observed. 
We believe it to be mostly due to the instability of 
isolated instantons with p.b.c. This conclusion is 
further supported by comparing the distributions 
(la,b) and (la',b') obtained from (la,b) by leav- 
ing out configurations of S' = Stt^ (see Fig. 4). 



(la'.b') sw=300 



3 



3 - 
2 - 
1 - 



0.2 



0.4 0.6 
p(fm) 



0.8 



Figure 4. Size distribution for (la) and (lb) once 
subtracted the configurations with S — Stt^. 



More detailed results, including larger statis- 
tics for /3 = 2.6 and an analysis for quantities 
like distance distributions, cross - correlations of 
size and distances and charge distributions will 
be presented elsewhere. 
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